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Preface 

A study to determine the feasibility of employing beamed electromagnetic energy for vehicle 
propulsion within and outside the Earth’s atmosphere was co-funded by NASA and the Defense 
Advanced Research Projects Agency that began in June 2010 and culminated in a Summary Presentation 
in April 2011. A detailed report entitled “Beamed-Energy Propulsion (BEP) Study” appeared in February 
2012 as NASA/TM — 2012-217014. Of the very many nuances of this subject that were addressed in this 
report, the effects of transferring the required high energy-density electromagnetic fields through the 
atmosphere were discussed. However, due to the limitations of the length of the report, only a summary of 
the results of the detailed analyses were able to be included. It is the intent of the present work to make 
available the complete analytical modeling work that was done for the BEP project with regard to 
electromagnetic wave propagation issues. In particular, the present technical memorandum contains two 
documents that were prepared in 201 1. The first one, entitled “Effects of Beaming Energy Through the 
Atmosphere” contains an overview of the analysis of the nonlinear problem inherent with the transfer of 
large amounts of energy through the atmosphere that gives rise to thermally-induced changes in the 
refractive index; application is then made to specific beamed propulsion scenarios. A brief portion of this 
report appeared as Appendix G of the 2012 Technical Memorandum. The second report, entitled “An 
Analytical Assessment of the Thermal Blooming Effects on the Propagation of Optical and Millimeter- 
Wave Focused Beam Waves For Power Beaming Applications” was written in October 2010 (not 
previously published), provides a more detailed treatment of the propagation problem and its effect on the 
overall characteristics of the beam such as its deflection as well as its radius. Comparisons are then made 
for power beaming using the disparate electromagnetic wavelengths of 1.06 pm and 2.0 mm. 
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Chapter 1. — Effects of Beaming Energy Through the Atmosphere 

1.1 Executive Summary 

Power beaming the large amounts of energy through the atmosphere that is required for propulsion 
either within the atmosphere or in transferring payloads from LEO to GEO will be met with the extremely 
deleterious effects of the thermal nonlinearities induced by atmospheric heating from beam power 
absoiption. The resulting phenomena, collectively called thermal blooming, will be the source of two 
major effects, i.e., beam steering away from the intended target vehicle as well as beam broadening which 
will make the beam larger at the target than intended. Both of these effects have been modeled and are 
shown here to be remedied by the use of high order (10 th order Elermite) phase compensation at the 
transmitter. Although it has been shown that these perturbations can be mitigated in principle, there may 
remain many engineering obstacles that must be overcome, such as elimination of mechanical jitter of the 
transmitter platform, etc., during beam propagation. Elowever, there is nothing in the prevailing physics of 
the situation that would preclude power beaming through the atmosphere as discussed in this report. 

It is important to begin to capture the prevailing effects and overall system operation using a scaled 
atmospheric experiment that would simulate the realistic environment in which an adaptive optics system 
must operate, from the variable wind velocity up to the mechanical jitter of the transmitter platform. It is 
recommended that, due to the immediate availability of high power sources, scaled atmospheric 
experiments be implemented (tailored after the Scaled Atmospheric Blooming Experiments (SABLE) 
Project by Lincoln Laboratory in the early 1990s) in which the operation of a closed-loop millimeter wave 
adaptive optics algorithm is assessed in the presence of a moving extended target. Additionally, this 
experimental scenario can also be used to address the issue related to the possibility of air ionization, and 
subsequent breakdown, across the apertures of the combined millimeter wave gyrotron sources as 
discussed elsewhere in this report. 

Finally, there is the need to look at various general beam wave profiles such as hypergaussian as well 
as fractional charge (in the topological sense) Laguerre-Gaussian beam waves that show great promise in 
their ability to be robust with respect to atmospheric nonlinearities. However, it is important that the 
modeling effort be kept to the level of yielding analytical results, rather than requiring numerical 
evaluation, so as to capture all the nuances of the physics involved and, at the same time, provide a tool 
for overall system evaluation as well as the design of adaptive optics algorithms. Finally, the model 
should be incorporated into a trajectory analysis program so a cadre of launch geometries can be 
evaluated from the point of view of atmospheric thermal nonlinearities. 

1.2 Introduction 

When delivering large amounts of power through the Earth’s atmosphere via millimeter or infrared 
‘beams’ (i.e., laser beams or beams formed at the output of a millimeter wave antenna system), many 
propagation mechanisms must be addressed that may be potentially deleterious such power transmission. 
The most obvious one is the ever-present random variation of the atmospheric refractive index due to 
local temperature variations known as ‘turbulence’. This naturally occurring phenomena is driven by 
thermal convection of heat from the Earth’s surface; once the resulting air motion exceeds a critical value 
of velocity, laminar flow essentially evolves into turbulent flow and fluctuations in the temperature 
distribution becomes statistically random (Ref. 1). These temperature fluctuations then act directly on the 
prevailing refractive index, thus rendering the refractive index a random quantity. These refractive index 
variations randomly focus and defocus the intervening electromagnetic wave field. Thus, the atmosphere 
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can be considered to be composed of ‘lenses’ of random focusing and defocusing characteristics that, due 
to the gross atmospheric motion due to wind, move across the beam. This gives rise to many beam quality 
variations; the major ones being beam broadening and beam steering. The statistical analysis and 
modeling of this type of atmospheric propagation as a long and rich history and has resulted in analytical 
descriptions for the impact of turbulence on the operation of systems relying on such beam propagation. 
Many models and descriptions exist for the ‘engineering analysis’ of the operation of transmission 
systems that rely on the propagation of electromagnetic beam propagation in the atmosphere (for a good 
recent treatment, see Ref. 2 and the references therein). 

The scenario discussed above may be considered as ‘passive’ electromagnetic wave propagation, i.e., 
the wave field moves through an atmosphere the refractive index of which is determined by other sources, 
not the field itself. However, as the energy density of the beam increases, absoiption of the beam energy 
by atmospheric gas components results in local heating of the atmospheric which does indeed act directly 
on the refractive index causing it to decrease in value. The possibility of this situation was first advanced 
in 1966 (Ref. 3). This thermal change of the refractive index field then acts on the electromagnetic wave 
field causing it to also change, and so on. The propagation scenario now becomes an ‘active’ one, 
whereby the propagating field modifies the very medium it which it exists. This heating process is called 
‘thermal blooming’ and substantially differs from that of the passive propagation discussed earlier 
(Refs. 4 and 5 (Ref. 5 contains a very comprehensive review of work and references that existed up to 
1990)). Here, a ‘thermal lens’ is created within the atmosphere by the heating due to the energy density of 
the beam. This ‘self-action’ of the beam will not only bend the beam into regions of higher refractive 
index (beam steering), but convection within the atmospheric fluid will also arise which is the source of 
self-induced turbulent flow of the medium. The situation is further complicated when one includes the 
effects of atmospheric wind and aerosols and the abovementioned passive propagation effects. 

Defocusing and other such associated nonlinear thermal blooming distortions of the beam cross-section 
will then result. In extreme cases of very large energy densities, the propagating beam will essentially 
break up into smaller beams, or filaments, which severely constrains the amount of energy density that the 
beam will be able to possess as it travels through the atmosphere. Unlike the situation of passive 
propagation, the thermal blooming mechanism introduces nonlinearities into the analysis of the 
phenomena that substantially complicates a complete mathematical description. Complete analyses of 
these types of propagation scenarios can only be done numerically, which was a major activity within the 
United States and Russia in the late 1980s. Other than the usual ‘order-of-magnitude’ estimates using the 
equations of fluid mechanics and wave propagation, only numerical modeling of the effects of 
atmospheric thermal nonlinearities abound in the literature. Analytical treatments appropriate for an 
engineering analysis of atmospheric propagation systems encountering thermal blooming have been 
lacking, especially those that endeavor to describe the result of adaptive correction of such nonlinear 
effects. This situation makes a comparative assessment of the operation of through-the-atmosphere power 
transmission difficult. 

In order to assess the atmospheric effects on the various beamed energy propulsion scenarios 
considered in this report, the propagation environment of high-energy electromagnetic wave transmission 
within the Earth’s atmosphere will be presented. First, a brief review will be given in Section 1.3 of the 
major physical mechanisms that prevail in the atmosphere that will deleteriously affect energy transfer via 
electromagnetic beam waves. The critical power thresholds for the wavelengths of 2.0 pm (infrared) and 
2.0 mm (millimeter wave) will be derived for various transmitter output aperture sizes at which these 
propagation mechanisms arise and must be addressed. Once this has been accomplished. Section 1.4 will 
present a propagation model that will be used to calculate the thermally-induced beam broadening and 
steering that will occur for a 2.0 pm laser beam, for a variety of aperture sizes and output powers, 
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transmitted from the Earth’s surface to a target at distances of 800 to 35200 km. The adaptive correction 
to mitigate the thermal nonlinearity effects on the beam will then be considered. Of the various 
performance parameters that an adaptive optics system can be designed to optimize at the target, it will be 
the minimization of the beam radius at target that will be the optimization parameter. Section 1.5 will then 
dwell on the dynamics that an adaptive optics system must satisfy, particularly with regard to the control 
delays inherent with the propagation distances that are involved. Finally, Section 1.6 will highlight some 
scaled ground experiments that can be performed using available high power millimeter wave sources. In 
particular, the effects of using an extended target with a closed-loop adaptive optics approach can be 
assessed. At this same time, the possible electrical breakdown of moist air in the vicinity of the combined 
gyrotron outputs can be studied. 

1.3 The Deleterious Atmospheric Propagation Mechanisms at High Energies 

At the power densities that are required for beamed energy propulsion, several aspects of propagation 
within the atmosphere must be addressed that can potentially perturb the beam wave and render it 
unreliable for energy transfer to a small target. Once these aspects have been identified, physically 
understood and mathematically modeled, appropriate mitigation procedures can be specified and applied 
to the problem to optimize energy transfer. These propagation features and the critical power levels at 
which they will be seen will now be given. 

1.3.1 Ionization and Electrical Breakdown 

When considering the transmission of high energy electromagnetic waves through the atmosphere, 
the first phenomenon that is considered is air ionization and the subsequent electrical breakdown (Ref. 6, 
a good introduction and the original work). Within the atmosphere, this occurs via the process of cascade 
ionization whereby a free electron is created by multiquantum absorption by the atmospheric gas and 
through the inverse bremsstrahlung process, accelerates and subsequently collides with an atom. The 
collision produces another electron, and both accelerate, collide with other atoms, and so on. This 
cascading process terminates in the release of light known as electrical breakdown. At a wavelength of 
K = 10.6 pm (i.e., the wavelength of high power CO 2 lasers), electrical breakdown occurs at a power 
density of « 10 9 W/cm 2 at atmospheric pressures and densities typical of those at sea level. This intensity 
is reduced by two orders of magnitude by the formation of the shock front due to the explosive detonation 
of atmospheric aerosols. Thus, for purposes of comparison, one can choose a power density of « 10 7 
W/cm 2 as the threshold for electrical breakdown within the atmosphere. 

At millimeter wavelengths in the 100 to 200 GHz frequency range, electrical breakdown occurs at a 
smaller power density of « 10 10 W/m 2 . One can then assume the worst case of a twofold decrease in this 
value due to atmospheric aerosols, thus giving a threshold power density of® 10 8 W/m 2 . 

1.3.2 Induced Molecular Polarization — The Kerr Effect 

The next process that must be considered is the Kerr effect, i.e., the process whereby the permittivity 
of the atmosphere, which is a function of the polarizability of the constituent gases, varies due the intense 
electric field of the laser beam acting to induce molecular orientation (Refs. 7 and 8). As discussed in the 
Appendix, the Kerr process can be described by a single differential equation, making its analysis 
straightforward as compared to the thermally-induced permittivity variations to be discussed in the next 
section. The associated relaxation times for molecular polarizability within the atmosphere are ~ 10 11 sec. 
so in the case of pulsed laser propagation, the medium can be considered in steady state for pulse lengths 
larger than this value. As the brief analysis in the Appendix shows, the threshold value P cr K of beam 
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power at which the Kerr effect is only a function of the Kerr constant and the wavelength (see Eq. ( A17)); 
it becomes an issue as a perturbing propagation mechanism at P cr K ~ 1 .5 x 1 0 9 W for the wavelength of 
X = 2.0 pm and P c ,k ~ 1.5><10 16 W at X = 2.0 mm. The Kerr effect serves to increase the permittivity in the 
location of the beam and thus results in what has become to be known as “self- focusing”, i.e., the beam 
tends to move into regions of larger permittivity and thus acts as if it has encountered a lens. The lens is 
essentially induced by the beam through the Kerr effect. This self- focusing phenomena results in a spatial 
instability within the beam causing a break-up of the beam into individual filaments, where each filament 
takes on the size which corresponds to P c ,k and the associated electric field strength. 

1.3.3 Induced Heating of the Atmosphere — Thermal Nonlinearities and “Thermal Blooming” 

The absorption of electromagnetic radiation by atmospheric gases becomes a source for the 
generation of heat. At large field intensities, the resulting temperature increase also changes the 
atmospheric permittivity. However, unlike the Kerr effect, the permittivity decreases and the beam tends 
to defocus (Refs. 4, 5, 9 to 12). Additionally, one must also admit the description of the atmospheric fluid 
dynamic processes that are elicited due to local heating and thus incorporate the Navier-Stokes and heat 
transfer equations for a moving medium; the medium is in motion not only due to the Archimedean forces 
that appear but also due to wind which serves to remove the heat that is generated. The problem in its 
entirety is a severely complicated one and is usually relegated to numerical analysis applied to specific 
cases. However, within certain approximations, analytical results can be obtained. As shown in the 
Appendix, one can obtain an expression for the critical power P c ,t the beam must possess above which 
thermally-induced propagation issues can become prominent. Due to the interplay of the several physical 
mechanisms that prevail in this propagation process, one must specify the prevailing wind velocity V as 
well as the effective radius r e ff of the output aperture of the beam source and the wavelength dependent 
absoiption constant a. As given by Equation (A25), one has for V= 4.4 m/s (~1 0 mi/hr), r e ff = 10 m and 
X = 2.0 pm, for which a ~ 4x 1CT 6 m _1 , P cr r = 16.7 W, a surprisingly small power. This is mainly due to 
the very large radius subtended by the beam through which the field-induced temperature increase must 
dissipate due to wind. For a more typical radius of r e ff= 0.1 m, one has P it t ~ 1.7 KW. At the much larger 
wavelength X = 2.0 mm, one has a ~ 2.2x1 0" 4 m 1 giving for r e /f= 10 m, P„t = 3.5xl0 5 W. 

In the case where the radiation source is pulsed with time duration of 1 ps, the threshold power 
becomes much higher due to the fact that the atmosphere does not have sufficient time to heat during the 
movement of the pulse through the atmosphere. As obtained from Equation ( A3 1), one has for the critical 
power at X = 2.0 pm, P cr p= 3.7xl0 7 W and for X = 2.0 mm, P cr p= 7.7x1 0 11 W. 

1.3.4 Atmospheric Aerosols 

The appearance of aerosols in the form of fog and clouds serve to be very efficient absorbers of 
electromagnetic radiation (Refs. 13 and 14). At the energy densities that prevail for beamed energy 
propulsion, aerosols will essentially explosively detonate thus diminishing their absoiption and scattering 
abilities. The dynamics of this complex process has been well studied and it has been established that, 
once a steady state has been established, propagation ‘channels’ will appear in the medium through which 
the beam can propagate (Refs. 15 and 16). 
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1.3.5 Atmospheric Turbulence 

The effect of random variations of the atmospheric permittivity field due to heating of the 
atmosphere, (the sources of which are solar absorption in the atmosphere and on the Earth’s surface) is 
independent of the beam power level for the lowest critical power thresholds derived above for the 
thermal blooming case. This ever-present deleterious mechanism has the most heritage in terms of study 
and understanding. The effects of this phenomenon on the propagation of a beam wave show themselves 
in terms of loss of spatial coherence, beam steering and beam broadening and are attendant with thermal 
blooming (Ref. 17). Of course, these characteristics are the same as seen with thermal blooming but the 
adaptive compensation in the case of turbulence is straightforward in that the effect being corrected is not 
the result of a nonlinear process. 

The results given in the above subsections are condensed in Table 1 and stated in terms of the critical 
power the appropriate propagation mechanism for output aperture sizes (where they matter) of 1, 10, and 
30 m. Also shown, for reference, are the general equations for the critical power and associated critical 
power density. As can be seen, discussions of the use of beam powers up to 1 00 MW at infrared 
wavelengths easily exclude the effects of electrical breakdown and the atmospheric Kerr effect leaving 
thermally induced nonlinearities and, of course, turbulence, to be dealt with. For millimeter wavelengths, 
the thresholds for electrical breakdown are three orders of magnitude smaller than for the infrared region. 
Hence, aperture sizes larger than 10m must be used for 1 00 MW power levels at these wavelengths using 
contiguous apertures. 


TABLE 1.— THRESHOLD POWER LEVELS BEYOND WHICH PROPAGATION EFFECTS OCCUR AT 
INFRARED AND MILLIMETER WAVELENGTHS FOR VARIOUS OUTPUT APERTURE RADII r e ff 


Propagation Mechanism 

X = 2.0 pm 

X = 2.0 mm 

General relation for critical power Per 
and power density p cr 

f 1 m 

Electrical Breakdown, r e g- = 1 10m 

[30m 

3.1xl0 n W 
3.1 x 10 13 W 
2.8 x 10 14 W 

3.1 x 10 8 W 
3.1 xlO 10 W 
2.8xlO u W 

No General Relation 

Kerr Effect 

1.5 x 10 9 W 

1.5 x 10 16 W 

aK 8>k 2 s K ’ PaK 8nk 2 r 2 ff s K 

Thennal Blooming (CW), 
f 1 m 
r eff = 10 m 
[30m 

Thennal Blooming (Pulsed, 1 ps) 

166.7 W 

16.7 W 
5.6W 

3.7xl0 7 W 

3.5xl0 6 W 
3.5xl0 5 W 
1.2 x 10 5 W 

7.7xlO n W 

n npC p V __ _ pC p V 

PcrT j 9 ? PcrT , 9 

k z r e ff£ T a k z r^z T a 

_ ttp C p _ pC p 

*crP ~ T 9 9 PcrP . 9 9 

k-ae T t p k z riv-as T t p 

Turbulence 

>0 

>0 
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1.4 Modeling the Effects of Atmospheric Thermal Nonlinearities on Beamed Energy 
Propulsion and the Required Level of Their Mitigation 

Although it is not the prevue of this work to physically model the thermal nonlinearities attendant 
with the propagation of high energy radiation through the Earth’s atmosphere, it has been found necessary 
to introduce a simple analysis employing the peculiarities of the propagation scenarios of beamed energy 
propulsion in order to analytically assess and compare the overall effects of the atmosphere in these 
various scenarios as well as evaluate the level of the corrective adaptive optics that will be required to 
make feasible the goals of beamed energy propulsion. Two such power beaming cases will be considered 
here, i.e., 1) orbit transfer from a low Earth orbit (LEO) to a geosynchronous Earth orbit (GEO) from a 
ground based laser transmitter operating at a wavelength of 2.0 pm in the infrared spectrum and 
2) beaming to a vehicle being launched within the Earth’s atmosphere using a ground based array of 
gyrotrons operating at the wavelength of 2.0 mm within the far infrared spectrum. As orbital transfer from 
LEO to GEO presents the worst-case situation in these two beaming cases due to the distances to be 
traversed, it will be used as the example to form an atmospheric propagation model upon which the 
overall operation can be assessed. It will then be shown how the model can be used to analyze the 
millimeter power beaming case for which the results will also be presented. 

The major property of the use of such propulsion to transfer payloads from LEO to GEO is, of course, 
the proportion of the diffraction length Ld associated with the beam wavelength and radius, to that of the 
Focal length F and the thickness H of the atmosphere. Here, one has that the diffraction length L c { = kr^j- 

(where k is the wave number of the radiation, r e ff is the effective radius of the transmitting aperture; see 
Appendix A) and F is such that Ld/F > H/Ld. This circumstance allows for the natural diffraction of the 
beam to be neglected within the region in which the beam is perturbed be induced nonlinear effects. This, 
along with the associated occurrence of the very large values of parameters describing the nonlinear 
interaction of the electromagnetic radiation with the fluid dynamics of the atmosphere (see the Appendix 
A), makes for the creation of an easily implementable model to describe the major beam parameters that 
are of interest, i.e., the radius of the beam at the target and the displacement of the beam from the target 
when no adaptive correction is applied. Once these have been secured, the corrections to the various 
orders of aberrations of the initial phase of the beam can be applied and their effect at the target be 
assessed. It is important to note here that only the effects of thermal nonlinearities of the atmosphere are 
considered. The additional effects of the ever-present turbulent fluctuations of the atmosphere will not be 
included here as this problem has much heritage and well understood and techniques for its mitigation are 
established. However, due to the very large power that the beams will possess that are considered here, it 
is expected that the induced convective velocity within the atmosphere in the beam path will swamp any 
turbulent velocity fluctuations. Thus, the turbulent mechanism affecting beam wave propagation at these 
power levels will be rendered negligible. The goal of the present analysis is to isolate the effects of the 
much more deleterious thermal nonlinearities on very long-range power beaming and as well as estimate 
the required level of adaptive compensation to make the orbital transfer mechanisms discussed elsewhere 
in this report realistic and feasible. 
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1.4.1 A Quick Overview of the Model and the Description of Beam Behavior at a Target in LEO 
and GEO due to Atmospheric Thermal Nonlinearities 

At the outset, the basis of the propagation model will be reviewed, leaving the details of the 
derivation to a forthcoming publication. The initial equations that will be employed are given by 
Equations (A22) and ( A26) of the Appendix A, viz., 

f)F" FT' i 

2 i + V \,E" + RyT'E" = 0, E" = E' exp(-aL d z'/2), = \E'\ (1) 

dz' p dx 


where the dimensionless nonlinearity parameter Rv is given by 


Rv = 


k 2 r eff £ T aP 0 

npC p V 


( 2 ) 


connecting the total incident power Po in the beam with the effective radius r e // of the transmitted beam 
and the prevailing wave number k = 2n/'h of the radiation of wavelength k. The atmospheric absoiption 
coefficient at the particular wavelength is ql,Zt=\ oz/oT\ is the variation of the permittivity with respect 
to the associated temperature variation T, p is the density of the atmosphere, C p is its specific heat at 
constant pressure, and V is the total relative velocity of the beam slewing across the atmosphere and that 
of atmospheric wind. ( All other variables and their normalization are introduced and discussed in the 
Appendix A.) Using the parameter values discussed in the Appendix A, one has, for a 1 MW beam, i.e., 

Po = l.OxlO 6 W, the value Ry= 6.0><10 4 W for X = 2.0 pm. In this instance, the second term in the first 
relation of Equation (1) can be neglected relative to the third term, i.e., the evolution of the beam within 
the atmosphere due to natural diffraction is negligible as compared to the phase perturbation due to the 
thermal nonlinearity acting on the permittivity function. Hence, the model equations become 

P)j7 n r)T f i 

2i— + R V T'E" = 0, E" = E'exp(-aL d z'/2), — = \ET (3) 

dz' dx 


The solution of these relations must now be augmented with the prevailing boundary condition for the 
initial radiation profile of the beam at the output aperture of the transmitter at z' = 0. Using a Gaussian 
beam profile, one has from Equation (All) 


£"(p',0) = Hoexp 




y' = \ + i^-, p '-=x' 2 +y' 2 
F 


(4) 


However, in most power beaming applications, the focal length F of the beam will be set to the distance L 
to the target, i.e., F = L. Additionally, the full normalized representation for a Gaussian beam wave must 
be employed in the last relation of Equation (3), viz, 


E f (p',z') = ^ exp 

l + iyz 




( Y ) 

n' 2 


v2 y 

(l + iy'r'J 

P 


( 5 ) 


Finally, one must consider the effective thickness H of the atmosphere. Although there are several 
models that can be used to represent the absoiption constant a with respect to the height within the 
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atmosphere (the most popular is an exponential variation with height), it proves expedient here to use a 
single value for H and a corresponding effective value for a. 

The solution of Equations (3) and (4) can now be considered in the case for LJF > FULd ; to be sure, 
one has if® 10 km, F < 35,200 km, and for the aperture size for 30 m, Ld ~ 2.8><10 9 m at infrared 
wavelengths. 

1.4.2 Application of the Foregoing to the Calculation of Beam Spread and Deflection due to 
Atmospheric Thermal Nonlinearities for Power Beaming to LEO and GEO 

Again, leaving the details concerning the derivations, etc., to a forthcoming NASA Technical 
Memorandum, one can employ Equations (3) to (5) to obtain expressions for the two very important 
‘performance parameters’ for power beaming to a target, i.e., the amount that the beam radius is widened 
due to the induced thermal de-focusing (beam spread) and the amount that the entire beam is steered away 
from the intended target due to the fact that the propagation path of the beam tends to favor regions of 
larger permittivity, i.e., steer away from regions of higher atmospheric temperature caused by absoiption 
(beam displacement). It is important to note that, as discussed in the Appendix A, the wind direction has 
been taken to occur along the x-axis of a coordinate system whose origin is situated at the transmitter 
aperture with beam propagation occurring along the z-axis. Also, the radius of the transmitter aperture at 
z = 0 is given by an effective radius r e j0) which is related to the corresponding beam waist radius Wo by 
r eff {0) = W 0 /j2. 

Because of the asymmetry introduced into the problem by the wind velocity V moving along the x- 
axis of the originally circular beam, there are two effective radii r i ,ff rX (L) and r e ffy(L), that characterize the 
beam in the x and y axes, respectively, at the distance L to the target vehicle. Again, this circumstance is 
brought about by the removal of heat from the beam channel along the x-axis. Elowever, heat is not 
convectively removed from the beam along they-axis; it diffuses in the y direction, which is a much 
slower process than convection. The result is that the de-focusing in they direction is much greater than 
in the x direction and the beam becomes elliptical as it leaves the atmosphere and continues to enlarge 
until it interacts with the target in space. 

The results of the calculation for the effective radii and the target are given by 

r eff,x {L) = R e ff + -® 2 N LV ) (6) 

r#A L ) = % +0.866 (/>») + 1.61 2®F (P„) (7) 


where 


®nlv(Fo) 


kztU.FlP^ 

4npC p Vr eff (0) 


( 8 ) 


is yet another nonlinear parameter that enters the problem and characterizes, within the thin phase screen 
approximation, the level of the thermal nonlinearity effect on the beam wave, of initial power Po, through 
an atmosphere of thickness H. (See Appendix A for the definitions of the remaining variables.) The 
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quantity R e jf(L) = L/(kr e fj( 0)) is what the radius of a beam focused on the target at a distance L would be if 
the atmosphere were not present, i.e., if transmission occurred entirely within a vacuum. (The reader is 
asked to excuse the rather excessive notations related to the various beam radii; the initial radius of the 
beam at the transmitter, denoted here as r e /f( 0), is usually stated in terms of the waist radius Wo that is a 

factor of V2 smaller than the former. However, it is desired here to show how the effective radii ‘evolve’ 
throughout the propagation process.) Thus, the quantities given by the radicals in Equations (7) and (8) 
are the factors by which the beam radii increase at a space home target due to the induced thermal 
nonlinearities within the atmosphere. 

Additionally, as discussed above, there occurs a deflection A (L) of the beam axis into the direction 
opposite to that of the wind. This is given by 


A(T) = 2R e ff (L)@ nlv (P 0 ) (9) 

The implications of these results to the behavior of a beam wave to targets at LEO (L = 800 km) as well 
as GEO ( L = 35,200 km) will now be given. The operating wavelength of the continuous wave case 
considered here is taken to be X = 2.0 pm and the atmospheric wind velocity along the x-axis is taken to 
be V x = 4.4 m/s « 10 mi/hr. Three transmitter aperture radii of r r .//(0) = 5 m, r e ff( 0) = 10 m, and /v//(0) = 30 
m are considered. Table 2 displays the beam radii that would occur at a target in LEO as well as GEO if 
the atmosphere were not present. 

Only a subset of these possibilities will be considered in what is to follow; to keep things realistic, 
r e ff{ 0) = 5 m at GEO and r<^( 0) = 30 m at LEO will not be subject to further analysis. 

Figure 1 to Figure 4 display the plots of r et f rX (L), r e ff- y (L) and A (L) versus Po for transmission to LEO 
and GEO for the various initial beam radii at the output aperture using Equations (6) to (9). In addition to 
the values quoted in Appendix A, the nominal value of H= 10 km has been used in conjunction with the 
average absorption coefficient a « 4x 10 6 m 1 . As can easily be seen in all the examples, the thermal 
action of the atmosphere severely impacts the integrity of the beam wave at the target locations. 


TABLE 2.— BEAM RADII AT LEO AND GEO TARGETS FOR VARIOUS 
OUTPUT APERTURE SIZES THE WAVELENGTH IS X = 2.0 pm 



r e /f( 0) = 5 m 

r e ff( 0) = 10 m 

r e ff( 0) = 30 m 

LEO (L = 800 km) 

0.1 m 

0.051 m 

0.017 m 

GEO (L = 35,200 km) 

4.48 m 

2.24 m 

0.75 m 


NASA/TM- 


-2015-218726 


9 



Beam Radii at Target vs. Power, L=800 km, 
r cf f(0)=5 m, V v =4.4 m/s 


; eff (L) (m) 

3.0 , 

2.5 

2.0 

1.5 

1.0 ; 

0.5 
0.0 ' 




(a) 


Induced Beam Deflection vs. Power, L=800 km, 
/- c ff(0)=5 m, V v =4.4 m/s 


Deflection (m) 



Figure 1. — (a) Beam radii (in m) along semi major and minor axes versus power (in MW) for a 5 m beam into LEO. 
(b) Beam deflection from target (in m) versus power (in MW) for a 5 m beam into LEO 
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Beam Radii at Target vs. Power, L=800 km, 
r ef f(0)= 10 m, V v =4.4 m/s 

r eff<L) < m ) 




(a) 


Induced Beam Deflection vs. Power, L=800 km, 
r eff (0)= 10 m, V X =4A m/s 

Deflection (m) 



TO (MW) 


Figure 2. — (a) Beam radii (in m) along semi major and minor axes versus power (in MW) for a 1 0 m beam into LEO. 
(b) Beam deflection from target (in m) versus power (in MW) for a 10 m beam into LEO. 
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Beam Radii at Target vs. Power, L=35200 km, 
r eff (0)=IOm,V ( =4.4m/s 



Induced Beam Deflection vs. Power, L=35200 km, 


r e ff(0)=10 m, V x =4.4 m/s 

Deflection (m) 



PO (MW) 

100 


Figure 3. — (a) Beam radii (in m) along semi major and minor axes versus power (in MW) for a 1 0 m beam into GEO. 
(b) Beam deflection from target (in m) versus power (in MW) for a 10 m beam into GEO. 
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Beam Radii at Target vs. Power, L=35200 km, 
r e ff(0)=30 m, V X =4A m/s 


r eff(L) (n>) 
30 r 

25 

20 

15 

10 

5 

0 


A = 2.0/i 



10 


20 


50 100 


TO (MW) 


(a) 


Induced Beam Deflection vs. Power, L- 3 5 200 km, 
r eff (0)=30 m, V x =4.4 m/s 

Deflection (m) 



P0 (MW) 


Figure 4. — (a) Beam radii (in m) along semi major and minor axes versus power (in MW) for a 30 m beam into GEO. 
(b) Beam deflection from target (in m) versus power (in MW) for a 30 m beam into GEO. 


Comparing the radii that result at the targets in LEO and GEO, due to a beam that has traversed 
thermal nonlinearities within the atmosphere, to those quoted in Table 2 shows the rather unsatisfactory 
but not unexpected results that deem beaming power to such locations as unfeasible without the proper 
compensation to the beam applied at the transmitter. The mitigation of the induced deflection of the beam 
can easily be removed by a tracker but the remaining aberrations that remain that result in the beam 
broadening can only be removed by appropriate phase compensation applied at the transmitter. 

1.4.3 Phase Compensation of the Effects of Atmospheric Thermal Nonlinearities 

The first research into the compensation of thermal blooming effects considered the application of a 
Zemike polynomial expansion of the transmitted phase front (Refs. 18 and 19). This technique provides 
for an optimal representation of the phase front of the beam wave so long as the cross-section of the beam 
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remains circular; after all, Zemike polynomials are, by design, orthogonal on a unit circle. The problem 
encountered in the above development is different; a beam asymmetry arises due to the differences in 
convective versus diffusion cooling of the atmospheric channel in which the beam propagates. The model 
briefly reviewed in Section 1.4.1 allows for the analytical assessment of phase compensation of the 
deleterious effects shown in Section 1.4.2. 

The different behavior of the beam along the transverse axes suggest the use of a more general set of 
orthogonal polynomials than those of Zemike. Here, since the initial form of the beam wave is given by a 
Gaussian function, the transmitted phase front S(x,y) at the output aperture of the transmitter can be 
expanded in a set of orthogonal polynomials associated with a Gaussian weight function, i.e., the Hermite 
polynomials H„(x) and H,„(y) (Ref. 20). Thus, 


M N 

S ( X ’T) = EEWf, (x)H m (y) 
n m 


( 10 ) 


The expansion coefficients a H are determined by applying one of several ‘performance metrics’. For 
example, for power beaming applications considered in this study, it is desired to shape the phase front so 
as to minimize the radius of the beam at the target so as not to have the beam interfere with the adjacent 
structures of the vehicle. Hence, the performance metric to be minimized through the selection of the a nm 
is related to the electric field E"(p',L) of the beam at the target vehicle by 


\\E"{p',Lfp'^dp' 

n (^ > P' = V( x ' 2 + T 2 ) dD 

]\E”(p',Lfp'dp' 

o 


in which E"( p',L) results from the application of the boundary condition 


£’"(p',0) = Hoexp 


v 2 , 




y f = 1 + i—^L 
L 


( 12 ) 


to the solution of Equation (3). Ideally, the number of terms of the expansion given in Equation (10) is 
infinite, i.e., N,M— > oo. The linear terms for which n + m = 1 determine the inclination correction of the 
beam wave; the quadratic terms for which n + m = 2 determine the focusing correction and the higher 
order terms for which n + m> 3 give the higher order aberration corrections. However, in realistic 
applications, the number of aberrations N and M is finite and, in fact, it is desired to find the smallest 
number of aberration corrections that can be used to represent the compensated phase. The specific 
expressions for the coefficients a„ m are found by equating to zero the derivative of H(Z.) with respect to 
these coefficients. 

The propagation problem defined here, as well as the associated minimization problem, can be solved 
analytically and will also appear in a forthcoming NASA Technical Memorandum. The results applied to 
the calculation of the phase-corrected beam radii are as follows. Due to the elliptical shape of the beam 
which is characterized by two radii at the target vehicle, it facilitates the calculations to consider a single 
effective radius r e V (T) defined by 
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Applying the modeling procedure as outlined above, one finds for the radius r' e ff Corr (L) of a focused 
beam wave corrected for the first M phase aberrations represented by the model of Equation (10) 


r ejf ,Corr (^) 7 ? eff (^) 


l + M»l LV (P„)-V->l LV (P 0 ) 


w 


yy /U0)M0) 

, „ 2 3(«+mV2„! m! 

7Z=I 772 =U 


^ + m)n\m\ 


71 

+ 4 ^ 

«=0V 




V' 


( 14 ) 


2 M2 n\ 


2 n n\n 


Figure 5 to Figure 8 show the results of this relation applied to the various FEO and GEO beaming 
cases considered in the last section. 


Beam Radii at Target vs. Power After Hermitc 
Polynomial Correction. L=8(K) km. 
r cff (0)=5 m, V / ,=4.4 m/s 

r cff.Corr< L > < m > 



0.0 P0 (MW) 

1 2 5 10 20 50 too 

Figure 5. — Beam radius (in m) at target versus power (in MW) for a 5 m beam into LEO for various orders of phase 

correction. 
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Beam Radii at Target vs. Power After Hermite 
Polynomial Correction, L=800 km. 



0.0 1 1 — ■■ ‘ - ' 1 — — 1 • ■ ■ P() (MW) 

12 5 10 20 50 100 

Figure 6. — Beam radius (in m) at target versus power (in MW) for a 10 m beam into LEO for various orders of phase 
correction. 


Beam Radii at Target vs. Power After Hermite 
Polynomial Correction, L=35200 km. 



Figure 7. — Beam radius (in m) at target versus power (in MW) for a 10 m beam into GEO for various orders of phase 
correction. 
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Beam Radii at Target vs. Power After Hermite 
Polynomial Correction. L=35200 km, 
r cff =30 m. V v =4.4 m/s 


r eff.Coir (L) <m) 
8 




Figure 8. — Beam radius (in m) at target versus power (in MW) for a 30 m beam into GEO for various orders of phase 
correction. 


These plots display the fact that for initial beam powers of 10 MW or less, phase correction using 
M= 10 order aberrations returns the radius of the beam to its desired value at the target vehicle. Beyond 
this power level, only the case shown in Figure 8 using a beam radius of 30 m at the transmitter can a 10 th 
order phase correction be effective up to about 70 MW. 

One assumption that has been prevailing in all these calculations is that the wind velocity Fis a 
constant along the entire propagation path. In reality, this is certainly not the case and, in fact, there may 
be regions in which the F« 0, i.e., a stagnation zone, and the convective heating scenario assumed here 
will transform into a diffusive heating one characterized by much lower critical power thresholds. This 
case will not be studied here, as it is the subject of a recent study (Ref. 21), but it must be kept in mind 
that a statistical description of the wind profile along the propagation path will be required in a more 
careful examination. 

Finally, there is the whole issue of the effect that the beam profile will have on the level of thermal 
nonlinearity perturbations. In fact, very high-energy laser beams are characterized by initial beam profiles 
that are ‘tubular’, i.e., have a local minimum of intensity in the center which increases toward the outer 
portions of the beam. Such beams, which can be modeled by hyper Gaussian profiles, tend to be 
characterized by higher critical power thresholds. This is due to the fact that the center of the beam, with a 
lower intensity, is thermally cooler than the periphery and thus tends to bend toward the center and 
remain stable. This effect should be studied for the particular high-energy sources that will be used as the 
design of the power-beaming scenarios evolve. 
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1.4.4 Using the Propagation Model for the Analysis of Millimeter Wave Beaming Within the 
Earth’s Atmosphere 

The methodology developed above for the infrared power beaming case for LEO/GEO orbital transfer 
can be applied to the millimeter wave case for power beaming to vehicles launched from within the 
Earth’s atmosphere. The wavelength of concern here is X = 2. Ox 10“ 3 m = 2.0 mm with an output aperture 
radius r e f = 50 m. This gives Ld = 7.9x 10 6 m. The propagation distances involved range from L- 20 km 
to L = 120 km. Although the beam will not be directly focused on the target vehicle (the dimensions being 
envisioned will require the beam to have a focal length just beyond the vehicle), for the purposes of this 
discussion, F « L. The effective distance within the atmosphere responsible for most of the absorption at 
these wavelengths is H& 1.8 km. Elence, the condition A//A> H/Ld is easily met. Additionally, there is a 
characteristic length L T that is a measure of the distance at which thermally induced diffraction effects 
will occur. This characteristic length can be defined using the length Ld as well as the nonlinear distortion 
parameter Ry given by Equation (A24). In fact, a detailed analysis similar to that given in Appendix A 
gives for the thermal diffraction distance 


l nr e//P C P V 
V s tP< 


(15) 


At the wavelength considered in this section, the associated absoiption coefficient is a « 2.2x10 4 nf 1 
which gives for the total transmitter power of Po = 800 MW (which will be considered below), 

Lt~ 73 km. Hence, thermal diffraction effects will certainly not prevail within the region H « 1.8 km or, 
for that matter, within the entire troposphere and the model constructed above can be adopted here. (This 
circumstance is solely due to the very large output aperture that is being used.) 

Equations (6) to (9) can now be used to consider four scenarios for a typical power-beaming launch: 

1) Horizontal beaming (0 = 0°) at L = 20 km 

2) Beaming at 0 = 30°, L- 50 km 

3) Beaming at 0 = 50°, L = 100 km 

4) Beaming at 0 = 90°, L = 120 km 


In each case, the target is to be illuminated in an area with a radius of R e fj{L ) ~ 1.8 m. Thus, a dynamic 
focus is to be employed; the focal length is determined by 



« =1± -L, *W> _, 


2L, 


Lr 


(16) 


where the sign is used in front of the radical for a beam focused behind the target. The results of 
applying the propagation model to these cases are displayed in Figure 9 to Figure 12. 
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Beam Radii at Target vs. Power, L=20 km, 0=0 deg, 
r cff =50 m, V,=4.4 m/s 

r cfr (L) (m) 

25 r / 


20 


15 


10 


0 

1(H) 


A = 2.0 mm 



ISO 200 


.MX) 


500 700 

(a) 


P0 (MW) 


Induced Beam Displacement vs. Power, L=20 km, 0=0 deg, 
r eff =50 m, V v =4.4 m/s 


Displacement (m) 



o- 

100 


150 


200 


300 

(b) 


500 


700 


P0 (MW) 


Figure 9. — (a) Beam radii (in m) along semi major and minor axes versus power (in MW) for a 50 m beam, (b) Beam 
deflection from target (in m) versus power (in MW) for a 50 m beam. 
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Beam Radii at Target vs. Power, L=50 km, 0=30 deg, 
r cff =50 m, V x =4.4 m/s 


r eff (L) (m) 



A = 2.0 mm 


2 

()t 1 ‘ — P0 (MW) 

100 150 200 300 500 700 



(a) 


Displacement (m) 


Induced Beam Displacement vs. Power, L=50 km, 0=30 deg, 
r cff =50 m, V x = 4.4 m/s 



OjO 1 • — ■ — ■ — • — 1 — — — • — PO (MW) 

100 150 200 300 500 700 

(b) 

Figure 10. — (a) Beam radii (in m) along semi major and minor axes versus power (in MW) for a 50 m beam, 
(b) Beam deflection from target (in m) versus power (in MW) for a 50 m beam. 
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Beam Radii at Target vs. Power, L= 100 km, #=50 deg. 


r ef f=50 m, V r =4.4 m/s 

r cff (L) (m) 



0 PO (MW) 

100 150 200 300 500 700 


Displacement (m) 
2.0 r 


(a) 

Induced Beam Displacement vs. Power, L=100 km. #=50 deg, 
r e ff=50 m, V x =4.4 m/s 



(b) 

Figure 1 1 . — (a) Beam radii (in m) along semi major and minor axes versus power (in MW) for a 50 m beam, 
(b) Beam deflection from target (in m) versus power (in MW) for a 50 m beam. 
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Beam Radii at Target vs. Power, L= 120 km, 0=90 deg, 
r e f f =50 m, V x =4.4 m/s 

r eff (L) (m) 




PO (MW) 


Displacement (m) 


Induced Beam Displacement vs. Power, L= 120 km, 0=90 deg, 
r eff =50 m, V x =4.4 m/s 



0.0 
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(b) 


500 
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Figure 12. — (a) Beam radii (in m) along semi major and minor axes versus power (in MW) for a 50 m beam, 
(b) Beam deflection from target (in m) versus power (in MW) for a 50 m beam. 


As can be seen from these results, the worst case is, of course, the horizontal beaming at an elevation 
angle 9 = 0° as displayed in Figure 9(a) and (b). Flere, the beam traverses 20 km of atmosphere. The beam 
severely broadens along they- axis, perpendicular to the wind velocity, since the heat that is generated is 
not convectively removed in this direction as it is along the x-axis. Related to this phenomenon is the fact 
that the beam is deflected in the x-direction up to 12 m from the target at a power level of 800 MW. Of 
course, as the elevation angle of the beam path increases, the distance within the atmosphere decreases up 
to the best case in which 0 = 90° as depicted in Figure 12(a) and (b). It must be cautioned, however, that 
the model used here, as mentioned earlier, does not include viscous and Archimedean forces that occur 
during the thermally induced motion of the atmosphere within the vertical column in which the beam 
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propagates. For situations in which large elevation angles are realized, these effects must be included in a 
future extension of the analysis. Finally, it must be noted that in order to follow the relatively close target 
vehicle, the beam will have a discernable slewing velocity as compared to that of the LEO/GEO transfer 
case discussed earlier. This slewing velocity may (depending on the relative wind velocity) tend to lessen 
the thermal effects. This will certainly be the case for the upward movement of the beam. As for the 
horizontal component of the slewing velocity, it can either detract or worsen the thermal effects. Such 
aspects of the problem can only be assessed via a simulation of the particular launch environment. 

Just as with the infrared beaming case, the adaptive correction of these deleterious effects is 
accomplished through the use of phase compensation at the transmitter output aperture. Again, because of 
the asymmetry induced within the structure of the beam, it is advantages to represent the corrected phase 
front as an expansion of Flermite polynomials. Figure 13 to Figure 16 show the results of applying 
various orders of Flermite correction to the millimeter wave beam. The goal is, of course, to return the 
beam radius to the desired value of 1.8 m. As expected, the 9 = 0° case requires the most compensation; 
correction up to 10 th order aberrations is needed, after the removal of the beam deflection via a tracker. 
The other cases only need correction up to 5 th order. 


Beam Radii at Target vs. Power After Hermite 
Polynomial Correction, L=20 km, 0=0 deg, 
r ef f=50 m, V x =4.4 m/s 

Corrected r c ff(L) (m) 


A = 2.0 mm 



2nd Order 
3rd Order 
8th Order 
I Oth Order 
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Figure 13. — Beam radius (in m) at target versus power (in MW) for a 50 m beam into GEO for various 
orders of phase correction. 


Beam Radii at Target vs. Power After Hermite 
Polynomial Correction, L=50 km, 0=30 deg, 
r e |f=50 m, V x =4.4 m/s 

Corrected r.i-pL) (m) 

4 

3 : A = 2.0mm 
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Figure 14. — Beam radius (in m) at target versus power (in MW) for a 5 m for various orders of phase 
correction. 
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Beam Radii at Target vs. Power After Hermite 
Polynomial Correction, L= 100 km. 0=50 deg, 

/• cff =50 m, V v =4.4 m/s 

Corrected (m) 
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A = 2.0mm 

2 ; 1 
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q[ 1 . 1 1 po (MW) 

100 150 200 300 500 700 

Figure 15. — Beam radius (in m) at target versus power (in MW) for a 5 m for various orders of phase 
correction. 

Beam Radii at Target vs. Power After Hermite 
Polynomial Correction, L= 120 km, 0=90 deg, 
r ef c=50 m, V x =4.4 m/s 

Corrected r ef | (L) (m) 
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Figure 16. — Beam radius (in m) at target versus power (in MW) for a 5 m for various orders of phase 
correction. 
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The additional nuance that this millimeter wave power-beaming scenario has over the LEO/GEO 
orbital transfer case is that the target vehicle will have a much larger apparent velocity in both the traverse 
and longitudinal directions than an object executing a LEO/GEO transfer. This circumstance brings to the 
fore the conditions of stability of the control algorithms which compose a closed-loop adaptive optics 
system. As mentioned above, the beam will also have a significant slewing velocity that, in some respects, 
eases the effects of thermal nonlinearities, but in others, places a burden on the operation of the adaptive 
optics control system. This will be discussed below in Section 1.5. 

It finally should be mentioned that for the power levels considered in this millimeter wave case, one 
should consider the mechanism of explosive detonation (Ref. 16) of small atmospheric aerosols that may 
traverse the beam, especially in the horizontal beaming case. This will give rise to a time dependent 
absorption coefficient that my indeed lessen the long-term heating effects within the atmosphere. 
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1.4.5 A Tool for the Quick Assessment of Nonlinear Thermal Effects on Power Beaming 

Appendix A introduced several expressions for threshold power levels at which thermal effects will 
become prevalent. Although these derivations were essentially based on normalization and dimensional 
analysis, they can be given credence by using the results obtained above. However, these critical power 
levels need to be related to a specific property or characteristic of a beam wave that will impact the 
performance of the power beaming system. As there are many such characteristics, there too will be as 
many corresponding critical power levels; there is no single universal quantity that can be stated to 
convey the specific impact thermal blooming will have on the entire performance of a beamed energy 
system. However, the most important characteristic in power-beaming situations is the variation of the 
beam radius due to the defocusing that occurs during the heating of the atmosphere. As shown by 
Equation (7) and, e.g., in Figure 1 to Figure 4, the radius of a beam in an atmosphere with convection is 
severely perturbed in the direction perpendicular to the wind direction. In order to attempt to secure an 
analytical estimate at which the radius of such a beam will begin to increase due to thermal blooming, one 
can expand Equation (7) and get 


r e(f,y it)- Reff(L) 1 + 


°- 866 L d 0 /p\, 

2 ^ u a£fI/oJ + 


(17) 


and require the value for Po at which 


— o/A' <°> eraP °" = 

° r NLVy ’ npCpVL 


2 L 


= 1 


(18) 


Setting L = H, one finally obtains for the critical power P crR at which the beam radius begins to expand 

npCpV 


PcrR- 10 


k 2 r ef f(0)e T a 


(19) 


which is just 10 P cr T as obtained in the Appendix A. One can now use this relation for millimeter 
wavelengths and compare to the infrared case. Figure 17 and Figure 18 show plots of Equation (19), as 
well as the associated power density P crR j nrP (0) , versus aperture radius for both cases. Here, the 

nominal wind velocity is taken to be V= 0.5 m/s. Of course, Equation (19) is simple enough to apply a 
statistical model describing the value for V. 

Evaluating Equation (19) using the numerical values of the parameters that do not depend on the 
wavelength, one has for this critical power 


PcrR- 200 


ar eff 


MW 


(20) 


where P c ,r is in MW and X and r e ff are in meters. (The units of a are m '.) These results clearly show that 
propagation at millimeter wavelengths within the atmosphere can sustain much larger power levels than at 
infrared wavelengths. 
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Critical Power and Power Density 
versus Aperture Radius 
for A=2.0 mm, V=0.5 m/s 
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Figure 17. — Critical power and power density at which beam radius is affected due to thermal defocusing versus 
initial aperture radius for millimeter wavelengths. 
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Figure 18. — Critical power and power density at which beam radius is affected due to thermal defocusing versus 
initial aperture radius for infrared wavelengths. 


1.5 Dynamics of Adaptive Phase Compensation for Power-Beaming Applications 

Of the two major approaches that can be used to implement phase compensation in order to mitigate 
atmospheric effects, viz., wave front ‘reversal’ or conjugation, and wave front control employing the 
optimization of a specific system metric for performance, the first presents a fundamental problem for the 
mitigation of thermal nonlinearities in that it relays upon a principle of reciprocity, a basic tenant of which 
is violated in thermal blooming applications. The use of such a principle essentially requires, through a 
wave front conjugation condition, that the amplitude and phase be corrected. However, it is much easier 
to correct the phase perturbations and not the associated amplitude perturbations. This is fine if amplitude 
variations are small, e.g., astronomical imaging. In thermal blooming scenarios, amplitude variations 
cannot be neglected. Using this approach in power beaming is not warranted. What is recommended here 
is the use of a method using the actual radiant intensity at the target. What used to be called aperture 
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optimization or ‘tagging’ (Ref. 22), has now come to be known as target-in-the-loop and its 
implementation is known as ‘gradient descent optimization (GDO) wave front control’ (Refs. 23 to 27) 
and is suggested for the power beaming applications addressed in this report. Here, the control rule is 
based on the direct optimization of an easily measured system performance metric, such as the radiant 
intensity at the target. 

The model governing the operation of such an adaptive optics system is given by first identifying a 
performance metric IT = n[S(5)] as defined by Equation (1 1) that is dependent upon the phase given by 

Equation (10) and the array a = \a nm ] of associated expansion coefficients (Ref. 24). The coefficients in 

Equation (10) were selected for modeling purposes based on the minimization of the beam radius on the 
target vehicle. Here, the optimal selection of the values of a„ m during the actual operation of the system is 
given by the control rule 


dO-nm (/) 

^ n 

" m dt 

in which x„ m are the time constants and y is the control (or ‘update’) coefficient. Here, J will be taken to 
be the size of the focal spot at the target vehicle. The intensity can be inferred by recording the reflected 
radiation form the target at a point separated from the transmitter so that a slightly different and 
unperturbed propagation path is used. It must now be established that such a control structure will 
dynamically operate in the two extreme cases of power beaming considered here, i.e., infrared beaming to 
a vehicle in GEO and millimeter wave beaming to a near Earth vehicle. 

Consider first power-beaming to GEO. For purposes of a brief analysis in which the required 
temporal characteristics of the adaptive optics system (and thus the temporal stability) are to be derived, 
one can consider the isolated case for which n + m = 2, i.e., according to the discussion in Section 1.4.3, 
the dynamic focusing of the beam wave. Thus, one has from Equation (10) for the associated (now time 
dependent) aberration coefficient 


dU(t) 


da„ 


( 21 ) 


a 2o(t)- 


m 


( 22 ) 


where, e.g., Hiix) ~ x 1 . For the case in which the focus is to be placed at the target at a distance L, i.e., 
F{t) = L, and the one has from Equation (21), 


a 2o(0 __ ^" + 


~~ + a 20 ( 0 ) 


exp 


(23) 


where t a = x 2 o L 2 d j [yLr j is the adaptation time of the adaptive optics system. This very important time 

constant is, although proportional to X20, inversely proportional to the distance to the target. This well- 
known property has it that the adaptive system will possess a faster response the farther the target from 
the transmitter. If this were all that the description of the adaptation system required, then the state of 
focusing on the target will monotonically approach the desired result LJL. However, there is a limiting 
factor placed on this circumstance by the delay inherent in overall system response do to the propagation 
time to and from GEO. When the propagation delay is allowed to enter into the control rule, Equation 
(21) now becomes 
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( 24 ) 


damn (0 _ dU(t-td) 
dt da nm 


where td is the delay time ( in this case, light time to and from the target). Expanding the right side of this 
equation in a Taylor series in td and keeping the first few terms yields 


8a nm (t) _ an(?) e 2 n(?) t 2 d e 2 n(t) 

dt du nm da nm dtd 2 du nm d tj 


(25) 


The corresponding characteristic equation q for the resulting differential relationship for aio ( t ) = exp(-qt) 
in the example considered above is given by 


S± = 


ta td — \ft% 2 t a td 7 


(26) 


shows that the system converges to a stable state so long as t a > td, i.e., the adaptation period of the system 
is longer than the delay time. For a target in LEO, td > 0.005 sec and for GEO, td ~ 0.24 sec (for the round 
trip time). If, on the other hand td > t a , the system does not converge to the stable state and 
lim a 20 (t)—> oo . Thus, the adaptive system must be continuously ‘tuned’ (by adjusting y in response to 

t — ><» 

changes in td as well as the evolving values for L in the case of a moving target) so that the condition 
t a > td is satisfied. 

To these considerations must also be added the temporal delay in the overall system response. This 
will not be discussed here as these considerations can be found in the literature. Suffice it to say that the 
delay that is met with in beaming to GEO can be easily factored into the operation of the adaptive system. 
It must be remembered that the distance to the target is always increasing in time for a beamed power 
scenario and this, above and beyond that of the propagation delay, must be properly treated in the design 
and specification of the adaptive optics control system. 

Work on the actual implementation of the close-loop control of the gradient descent algorithm of 
Equation (21) has rapidly progressed. The technology is now mature and form s the basis for reliable 
adaptive control. The only possible nuance to its application is for power beaming to locations within the 
Earth’s atmosphere for surface or near-surface launching. Here, the target vehicle, unlike the LEO/GEO 
transfer case, is rapidly moving and is considered an extended target, i.e., its spatial extent cannot be 
relegated to a point object and its possible rotation and orientation becomes an issue for use as a reflector 
of radiation bask to the adaptive optics receiver. Here, the randomly rough surface of the target reflector 
creates a speckle field characterized by bright and dark regions of intensity has been found to complicate 
the implementation of a closed loop adaptive optics scenario. This can lead to the slowing of the 
convergence speed of the system and, ultimately, to the efficiency of its operation in terms of the 
optimization of the focal spot at the target. Methods using the pre-compensation of the beam show 
promise (Ref. 28). 

Although it has been shown above that the basic physics of the situation does not preclude beaming 
energy to a prescribed area on a target vehicle in GEO, the engineering implementation will offer some 
major obstacles to be negotiated. In particular, the beaming into a GEO location will require 30 
nanoradian pointing accuracy on the part of the transmitter optics. The presence of mechanical jitter of the 
transmitter platform must be continuously mitigated. Additionally, on the other side of the spectrum, 
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launching from within the atmosphere will challenge the application of the adaptive optics at millimeter 
wavelengths. The fact that the beam is rapidly slewing due to the relatively rapidly moving object makes 
it difficult for the closed-loop adaptive optics system to converge to a stationary value since the medium 
within the column of atmosphere in which the millimeter wave beam exists is constantly being 
exchanged. Here, instead of employing the target as the beacon source, the use of an artificial beacon 
placed ahead of the moving vehicle may help in characterizing the atmosphere ahead of the beam. This 
will assess the phase perturbations of the nonheated atmosphere but the problem still remains concerning 
just how the atmosphere will respond to the heating from the beam as it arrives at that particular column. 
This suggests that, due to the availability of high power sources at millimeter wavelengths, one configure 
an in situ experiment to assess the operation of both closed-loop versus artificial beacon adaptive optics 
approaches. This will be discussed further below. 

1.6 Experiments Assessing the Compensation of Thermal Nonlinearities of High Power 
Beams in the Atmosphere 

Experiments dealing with the induced effects due to thermal nonlinearities elicited by high power 
propagation are usually performed in a laboratory using liquids placed in cells in which the thermal 
nonlinearity thresholds are much smaller than in air. In these scaled laboratory experiments, where 
atmospheric turbulence effects are simulated by transmission phase screens, various adaptive optics 
algorithms have been tested and evaluated (Refs. 29 and 30). Some experiments have been performed in 
the open atmosphere along horizontal paths in a program called the Scaled Atmospheric Blooming 
Experiments (SABLE) directed by Lincoln Laboratory (Ref. 31). Work is also continuing along these 
lines in other countries (Refs. 32). All these experiments endeavor to evaluate the in situ operation of 
adaptive optics systems on a scaled basis. The same must be recommended for the operation of a beamed 
power scenarios discussed above. Due to the maturity and availability of high power sources in the 
millimeter range, it is recommenced that a scaled atmospheric experiment be performed on a moving 
target to assess the operation of various adaptive optics algorithms. Such an experiment should be 
modeled after the SABLE project. It is important to note that there currently does not exist a database that 
addresses the power beaming cases considered here. It is important to begin to capture the prevailing 
effects and system operation using a scaled atmospheric experiment which would simulate the realistic 
environment in which an adaptive optics system must operate, from the variable wind velocity up to the 
mechanical jitter of the transmitter platform. Additionally, this experimental scenario can also be used to 
address the issue related to the possibility of air ionization, and subsequent breakdown, across the 
apertures of the combined millimeter wave gyrotron sources as discussed elsewhere in this report. 

1.7 Conclusions and Recommendations 

The various deleterious propagation mechanisms associated with high-power electromagnetic wave 
propagation through the atmosphere have been discussed. In addition to turbulence, thermal nonlinearities 
associated with the absoiption of radiation by the atmospheric gas will contribute to the major effects of 
beam wave propagation for beamed energy propulsion at the power levels considered here. The simplified 
propagation model advanced here shows that the beam radius as well as its deflection are severely 
affected by the phenomena of thermal blooming. However, the model also showed that appropriate phase 
compensation at the transmitter output aperture can mitigate these effects and return the propagation 
situation to one that is acceptable for the power transfer requirements that are needed to be satisfied for 
beamed propulsion. For infrared transmission out of the atmosphere for LEO/GEO beaming, up to 10 th 
order (once the tilt has been removed) aberration correction will be needed in order to maintain a minimal 
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focal spot at a GEO location using a 30 m diameter transmitter aperture. For millimeter wave beaming 
within the atmosphere, 5 th order aberration correction will suffice with the exception of the horizontal 
case in which, once again, a 10 th order correction would be requires. In principle, these corrections allow 
the beaming system to operate within the prevailing specifications. In practice, however, some challenges 
prevail in the implementation. 

The adaptive optics approach recommended here for the LEO/GEO launch case is a closed-loop 
system that uses the specular reflection from the target vehicle as a ‘beacon’ source (i.e., a target-in-the- 
loop system). The target vehicle is seen by the adaptive optics system as a simple point reflector. The 
round-trip propagation time delay inherent in this scenario can be tolerated so long as the adaptation time 
of the adaptive optics system is properly set to a slightly larger time than that of the delay to assure proper 
convergence. 

In the case of millimeter wave beaming to a moving vehicle within the Earth’s atmosphere, the object 
is now considered as an extended target the rotation and orientation of which can complicate application 
of the closed loop system. The reflection from the vehicle will posses a speckle structure that can severely 
impact the wave front sensor used by the adaptive optics system. Here, it may be that the simpler artificial 
beacon method can be used but there is still an additional complication, viz., the rapid movement of the 
beam across the atmosphere. The changes induced in the beam column through the atmosphere will not 
all be those due to the adjustment by the adaptive optics and the system will not be able to properly adapt. 
This situation presents itself for a scaled atmospheric high-power millimeter wave transmission 
experiment in which both closed loop and artificial beacon-based adaptive optics systems are tested and 
evaluated. 

The propagation modeling for the various beamed propulsion scenarios presented in this report should 
be extended beyond that employed here in Section 1.4 by incorporating additional fluid mechanical 
descriptions of the atmosphere as well as more general beam wave profiles such as hypergaussian as well 
as fractional charge (in the topological sense) Laguerre-Gaussian beam waves which show great promise 
in their ability to be robust with respect to atmospheric nonlinearities (Ref. 33). However, it is important 
that the modeling effort be kept to the level of yielding analytical results, rather than requiring numerical 
evaluation, so as to capture all the nuances of the physics involved and, at the same time, provide a tool 
for overall system evaluation as well as the design of adaptive optics algorithms. Finally, the model 
should be incorporated into a trajectory analysis program so a cadre of launch geometries can be 
evaluated from the point of view of atmospheric thermal nonlinearities. 
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Appendix A. — The Derivation of the Fundamental Nonlinear Propagation 
Equations for High Energy Transmission Through the 
Atmosphere and the Associated Critical Powers 

For the propagation of scalar electromagnetic waves through a medium characterized by a variable 
permittivity As, one has for the equation describing the evolution of the electric field strength 


f 


2 ik 


dE_ 

dz v dt 


+ V 2 E + k 2 
P 


( As^ 


V ^0 


E = 0 


(Al) 


where p is a position vector in the plane transverse to the direction of propagation along the z axis, 
k = In IX is the wave number of the radiation field of wavelength /,, u is the velocity of the wave field in 
the medium ( u = c where c is the velocity of light), and so is the nominal permittivity (for the 

atmosphere, So « 1). As will be discussed below, the function describing As is, in general, dependent on 
the value of E, i.e., As = As (E), thus making the propagation problem a nonlinear one; only in the case of 
weak fields (the typical situation for atmospheric optics) is As independent of E. To simplify the analysis, 
one can, to a very good approximation, assume that u > /,,/x where / is the longitudinal size of the spatial 

variation of As and x is its temporal variation; this will preclude the description of the propagation pulses 
of very short temporal duration where dispersion will become important but will suffice for the purposes 
of this presentation. Thus, Equation (Al) reduces to 


2 ik—+viE+k 2 
dz p 


As 


E = 0 


l e o ) 


(A2) 


which is known in the literature as the quasi-optical approximation. The first two terms of Equation (A2) 
describe the propagation (within the quasi-optical approximation) of a beam wave (see Eq. ( A4)) in a 
vacuum; the last term is the perturbation added to the evolution of the beam wave through a medium 
described by the factor As. It is through this factor that the various effects of the atmosphere on 
propagation come into play. 

For the power ranges of interest in high energy laser propagation, the total permittivity function can 
be written as 


As(|e| 2 ,tJ = /So J + As k (|e| 2 ) + s t T (A3) 

in which a is the wavelength dependent absoiption coefficient; absoiption within the propagation medium 
gives rise to the increase in the associated temperature T from the nominal value, which, in turn, modifies 
the value of the permittivity by the amount stT where Sr = ds/dT is the variation of the permittivity with 
respect to T. (Strictly speaking, the temperature increase T is also a function of \E \ 2 but, as will be 
discussed below, it is not a direct function as is the Kerr effect.) This gives rise to thermal nonlinearities 
and the associated thermal ‘blooming’. Also, due to the large value of E in cases of very high energy 
propagation, induced orientation of anisotropically polarized molecules of the atmosphere occurs due to 
the interaction of the induced dipole moment of the molecule due to the intense electric field. Thi s 
contributes another perturbation Asa' to the permittivity through the quantity \E\ 2 . This phenomenon is 
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known as the Kerr effect. The functional relations of both Ae K ||/T| 2 j and T to the quantity \E\ 2 are also 

given through the prescription of a differential equation. Each of these contributions will be separately 
discussed below and the critical field powers at which they will occur within the atmosphere will be 
derived. Before this is done, however, it will prove to be convenient to make the following digression 
concerning the description of a beam wave. 

The consideration of the general solution for Equation (A2) makes it necessary to have an expression 
for the associated boundary condition. This brings to the fore the normalization of spatial coordinates. 
The boundary condition is given by the expression for the general form of a Gaussian beam wave, viz., 


E{p,z) 
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f ky' 

| P 2 


{ 2 , 

/ 1 + iyz 


where A o is the initial amplitude of the field and 


y = yi+/y2 ’ Yi =^f y2 4 


(A4) 


(A5) 


with Wo being the waist radius of the beam at the exit of the output aperture and F is the focal length. The 
initial condition for the field at the output aperture is thus given by 


£(p,0) = Aq exp 
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At this point, it is convenient to normalize the spatial coordinates with respect to the parameters of the 
problem. For example, if one defines the ‘effective’ radius of the beam wave at the exit of the output 
aperture as r e jj- = W t) j \[2 , (this definition is consistent with the definition for the power distribution in 

the transverse plane) one can define the dimensionless transverse coordinate p' = p/'/V// and write 
Equation (A6) as 


£(p',0) = zloexp 


v2 y 


p ' 2 + i 


r kW 2 


■.'2 


(A7) 


One can now apply the same coordinate normalization to Equation (A2). That is, since 




(AS) 


Equation (A2) can be written (letting so = 1 for atmospheric scenarios) 


2 dz p 2 


(A9) 
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Now, normalizing the longitudinal coordinate z with respect to the quantity defined by L d = k Wq jl, i.e., 
z' = z/Ld, and noting that 


dE(p,z) _ 2 dE(p',z') 

dz kWfi dz' 


Equations (A7) and (A9) become 


E (p',0) = Aq exp 



(All) 


and 

2/|^ + Vj.E + k 2 r 2 f As(|A| 2 ,T^E = 0 (A12) 

This form of the quasi-optical equation will facilitate the developments that follow. The quantity Ld is 
called the diffraction length and, as seen from Equation (All), is the length that the focus must achieve 
for the associated diffraction to become appreciable. 


A.l The Critical Power for the Atmospheric the Kerr Effect 

To isolate the effect of the Kerr nonlinearity from the thermal effects due to absoiption, one simply 
lets As = Asa (|A| 2 ) in Equation (A 12) and obtains 

2/0 + V?A + k 2 r 2 f Az K (|if ) E = 0 (A13) 

The differential equation that governs the evolution of Asa is given by 

T^A + A ej ,=s*|£| 2 (A 14) 

where t is the relaxation constant for the Kerr process and Sa is the Kerr constant for the atmosphere; 
typically, x ~ 10" 11 sec. Relative to the other processes to be considered (and the fact that very short pulse 
lengths are not being considered here), one can ignore the first term on the left of Equation (A14). In this 
case, Asa = saIAI 2 and Equation (A13) becomes 

f)F o 

2i— + Wl,E + k 2 r^ K \E\ 2 E = Q (A15) 

In the literature. Equation (A15) is sometimes stated to describe a ‘cubic’ medium due to the occurrence 
of E in the last term. Connecting \E\ 2 to the associated intensity / of the beam using the relation 
| E\ 2 = ( 8n/c)I and then to the corresponding beam power Po using Pq = nrk / , Equation (A15) becomes 
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(A 16) 


2 i — + V?, E + R k E = 0 
dz' p 

where Rk = PJPcrK is the nonlinear parameter (also called the distortion parameter) for the Kerr 
Nonlinearity defined in terms of the critical power P c ,k associated with the Kerr nonlinearity, 


PcrK ~ 


C 

Sk 2 sx 


(A17) 


For Earth’s atmosphere at sea level, ek~ 2.5x10 16 CGSE units (Ref. 8). Thus, for a laser wavelength of 
X = 2.0 pm, one has that the Kerr effect has a threshold at P c ,-k~ 1.5x 10 9 W. This value is less than one 
order of magnitude from the powers that are being considered (~70xl0 6 W, giving Ii K ~ 0.05) for power 
beaming applications. Hence, the nonlinearities associated with the Kerr effect may become a minor 
factor in power beaming applications (and can be easily removed via adaptive optics) but, as will be 
shown below, Rk is several orders of magnitude below the corresponding parameter for thermal 
nonlinearities which dominate the beaming process in the atmosphere. 

In the case of millimeter wave propagation at X = 2.0 mm, one obtains P c ,-k = 1.5xl0 16 W and at 
Po = 70x1 0 6 W, Rk ~ 4.6x1 0“ 9 thus alleviating any concerns about the Kerr effect at these wavelengths. 


A.2 The Critical Power for Atmospheric Thermal Effects 
A.2.1 Continuous Wave Source 

Using Equation (A3) in Equation (A 12) and neglecting the Kerr contribution gives 


2i— + S7lE + k 2 rl f 
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The temperature variation from the nominal atmospheric temperature is given by the application of the 
conservation of energy to the atmospheric heat budget, given in its entirety by 


pC p 


— + V-VT 

dt j 


- k v 2 r + a / 


(A 19) 


where p is the density of the atmosphere, C p is its specific heat at constant pressure, k is its thermal 
conductivity, V is the total relative velocity of the motion (slewing) of the beam across is the atmosphere 
and that of the atmospheric wind, and a is the absorption coefficient for the radiation of intensity 
I = (c/Stt) | E\ 2 . (Strictly speaking, one should also include the set of Navier- Stokes equations in the 
Boussinesq approximation to account for the viscous and Archimedean forces that occur during the 
thermally induced motion of the atmosphere. However, for the purposes of this discussion and subsequent 
analysis, a description at this level is not required.) Equation (A19) can be simplified straight away by 

first specializing the slewing and wind velocity to be along the x axis of the coordinate system, V = Vx , 
where x is the unit vector. Further, the diffusivity V 2 T will be taken to have contributions only in the 
direction transverse to the laser beam propagation, i.e., V 2 T « V?T . Finally, only the stationary heating 

case will be considered for this particular discussion whereby the time derivative can be dropped. (That is, 
the CW radiation is taken to be acting long enough for the steady state case to be achieved. This is not the 
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case for pulsed radiation which will be discussed below.) Applying these considerations to Equation 
(A 19) allows it to be written 


K-= x v r + 


a c 

8npC p 



(A20) 


where y = K/p C p is the associated thermal diffusivity. As was done earlier with Equation (A2), one can 
normalize the transverse coordinates with respect to r e /f. Additionally, the field value E can be normalized 

with respect to E' = E/Eo where Eq = yj8nl/c - ^J%Po ~J r ejf c (this is what was essentially done to 

Equation (A15) to obtain Equation (A 16)). Equation (A20) then becomes 
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Now, the two heat transfer mechanisms that are described here can be isolated. Consider the ratio of the 
coefficients of the two terms within the brackets of Equation (A21), i.e., q = Vr e gly. This dimensionless 
ratio is known in fluid mechanics as the Peclet number. Taking the nominal wind velocity V= 4.4 m/s 
(10 mi/hr), r e /f= 10 m and using the documented value for the thermal diffusivity of the atmosphere 
% = 2.12x10 5 m 2 /s, one has that iq = 2.1x] 0 fi > 1 allowing one to neglect the second term in the brackets 
In this instance, one can write 


n P C p r e ff V 

aP n 


dT 
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— = |E'| 2 
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(A22) 


where the temperature has been normalized with respect to the characteristic temperature 
Tv= aPo/npCpr e ffV. Finally, applying this prescription to Equation (A 18) gives 

8F' 

2 i + yl,E' + iL d aE' + k 2 r} ff z T T v T'E’ = 0 (A23) 

dz p 


As done with Equation (A 15), one can define the distortion parameter for the thermal nonlinearity 


R v = k 2 r 2 f e T T v = 


k 2 r eff e T aP 0 

npC p V 
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(A24) 


where 


PcrT = 


npCpV 
k 2 r e jj s T a 


(A25) 


is the critical power associated with the thermal nonlinearity in the case of large Peclet numbers. Using 
Equation (A24) and noting that the third term in Equation (A23) can be transformed away, one has 
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E " = E' exp(-aL rf z'/2) 


(A26) 


2/ — + VIE" + R V TE" = 0, 
dz' v 


The problem of propagation through a thermal nonlinearity thus reduces, to within the 
approximations adopted, to the solution of Equations (A22) and (A26). For the Earth’s atmosphere, 
p = 1.2 kg/m 3 , Cp = l.OxlO 3 J/K-kg, and &r~ -2.3x10 6 K _1 . (That is, the atmospheric permittivity 
decreases with an increase in temperature. This indicates that the laser beam will steer toward regions 
with larger permittivity, i.e., cooler regions. In the calculations, the absolute value of this quantity is 
used.) At X = 2.0 pm, a ~ 4xl0 _6 m _1 and using the earlier stated values for the other quantities, one finds 
that P crT = 16.7 W. Thus, for Po = 70x1 0 6 W ,Rv~ 4.2x1 0 6 . The very small value for the critical power 
occurs due to the rather large beam radius. The temperature must diffuse over the large cross section of a 
cylindrical column. Of course, larger values of V will increase P IT r. For a smaller, more typical value for a 
beam radius, r e /f= 0.1 m, one finds that P c ,t= 1667 W and Rv~ 4.2x1 0 4 for Po = 70 MW, still a rather 
large value for the nonlinearity parameter. As a comparison, a more typical laser power of Po= 100 KW 
corresponds to Ry « 60. It is thus seen that nonlinearities due to thermal effects at infrared wavelengths 
within the atmosphere will dominate those of the Kerr effect. 

The small values for P cr r at the large output aperture radii and the corresponding very large values for 
Rv do indeed have severe ramifications for the long distance propagation of beam waves at Po ~ 70 MW 
as shown earlier. However, these values facilitate the use of approximation procedures to be used in 
solving Equations ( A22) and ( A26) from which specifications for the adaptive correction of the various 
aberrations are derived and it is shown that complete mitigation is possible. 

The situation for millimeter wave propagation is better. For X = 2.0 mm, a = 2.2x1 0‘ 4 m' 1 and one has 
P cr T= 3.5x10 s W and for Po = 70 MW, AT = 200.8. Although the absorption coefficient a is two orders of 
magnitude larger than that for infrared wavelengths, the overall optical effects due to the thermally 
induced refractive index variations are much smaller due to the larger wavelength; the ‘thermal lens’ that 
is induced has weaker focusing capabilities at millimeter wavelengths. 

A.2.2 Pulsed Source 

In the case of a source that produces radiation pulses of temporal duration t p , the prevailing fluid 
dynamics of the situation differs from that of a CW source. Here, for t p < r e ff! V, the convection/diffusion 
equation, Equation (A19) is approximated as 


pC p — = al (A27) 

Coupling this equation with the propagation equation, Equation (A18), and following the same 
normalization procedure as outlined above for the CW case and normalizing the time to t p , one obtains in 
place of Equation (A26), 


c)F" 

2 i — + VlE" + RpTE" = 0, 
dz' p p 


E" = F'exp(-aP t / z'/ 2) 


(A28) 


NASA/TM— 20 15-218726 


36 



where 


^P c P r eff ^\dT _8T 

aPotp dt' dt' 


and the prevailing nonlinear parameter is now 


_k 2 & T at p P 0 

Kp = 

ftpCp 


Jb_ 

PcrP 


in which the corresponding critical power is 


^crP ~ 


^P C p 
k 2 as ftp 


(A29) 


(A30) 


(A31) 


which is independent of r e //. 

Taking a nominal pulse width of t p « 1 ps, one has for X = 2.0 pm, P cr p = 3.7xl0 7 W and for 
X = 2.0 mm, P cr p = 7.7xlO n W. Thus, temporally modulating the source does not elucidate the convective 
thermal effects that are inherent in the CW case and allows for much larger critical values of peak power. 
Pulse lengths smaller than 1 ps require the use of Equation (Al) which takes into account pulse dispersion 
through the medium. This effect is especially important when atmospheric turbulence is taken into 
account. However, for the purposes of the present analysis, which is to obtain the critical powers at which 
various nonlinear propagation effects show themselves, these effects are not considered. 

The issue of the repetition rate of a pulse train within the atmosphere in conjunction with scanning 
(slewing) the beam through a region of the atmosphere to aid in heat dispersion (for high pulse repetition 
rates) is characterized by yet another critical power threshold which will be larger than that given by 
Equation (A31). 

A.3 Conclusion 

The parameters Rk, Ry and Rp provide for a comparative assessment of the nonlinear propagation 
mechanisms that are induced due to the large field power densities. Thus, at the power levels considered 
here, the effect of Kerr nonlinearities are almost negligible relative to that of thermal effects and it is the 
later that must be completely understood so as to be effectively mitigated via properly designed adaptive 
optics algorithms. They do not, however, give an assessment of the behavior of specific laser beam 
properties, such as beam steering and beam broadening that accompany propagation. This can only be 
done through modeling using the specific propagation geometries that are being considered. As 
mentioned above, the very large values for R v in the case of the power levels considered for beamed 
propulsion are a benefit in that they allow for the use of approximation procedures to be applied to 
analytically evaluate the propagation and adaptive optics models. 
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Chapter 2. — An Analytical Assessment of the Thermal Blooming Effects on 
the Propagation of Optical And Millimeter- Wave Focused Beam 
Waves For Power Beaming Applications 


2.1 Introduction 

When delivering large amounts of power through the Earth’s atmosphere via millimeter or infrared 
‘beams’ (i.e., laser beams or beams formed at the output of a millimeter wave antenna system), many 
propagation mechanisms must be addressed that can be potentially deleterious such power transmission. 
The most obvious one is the ever-present random variation of the atmospheric refractive index due to 
local temperature variations known as ‘turbulence’. This naturally occurring phenomena is driven by 
thermal convection of heat from the Earth’s surface; once the resulting air motion exceeds a critical value 
of velocity, laminar flow essentially evolves into turbulent flow and fluctuations in the temperature 
distribution becomes statistically random (Ref. 1). These temperature fluctuations then act directly on the 
prevailing refractive index, thus rendering the refractive index a random quantity. These refractive index 
variations randomly focus and defocus the intervening electromagnetic wave field. Thus, the atmosphere 
can be considered to be composed of Tenses’ of random focusing and defocusing characteristics that, due 
to the gross atmospheric motion due to wind, move across the beam. This gives rise to many beam quality 
variations; the major ones being beam broadening and beam steering. The statistical analysis and 
modeling of this type of atmospheric propagation as a long and rich history and has resulted in analytical 
descriptions for the impact of turbulence on the operation of systems relying on such beam propagation. 
Many models and descriptions exist for the ‘engineering analysis’ of the operation of transmission 
systems that rely on the propagation of electromagnetic beam propagation in the atmosphere (for a good 
recent treatment, see Ref. 2 and the reference therein). 

The scenario discussed above may be considered as ‘passive’ electromagnetic wave propagation, i.e., 
the wave field moves through an atmosphere the refractive index of which is determined by other sources, 
not the field itself. However, as the energy density of the beam increases, absoiption of the beam energy 
by atmospheric gas components results in local heating of the atmospheric which does indeed act directly 
on the refractive index causing it to decrease in value. This thermal change of the refractive index field 
then acts on the electromagnetic wave field causing it to also change, and so on. The propagation scenario 
now becomes an ‘active’ one, whereby the propagating field modifies the very medium it which it exists. 
This heating process is called ‘thermal blooming’ and substantially differs from that of the passive 
propagation discussed earlier. Here, a ‘thermal lens’ is created within the atmosphere by the heating due 
to the energy density of the beam. This ‘self-action’ of the beam will not only bend the beam into regions 
of higher refractive index (beam steering), but convection within the atmospheric fluid will also arise 
which is the source of self-induced turbulent flow of the medium. The situation is further complicated 
when one includes the effects of atmospheric wind and the abovementioned passive propagation effects. 
Defocusing and other such associated nonlinear thermal blooming distortions of the beam cross-section 
will then result. In extreme cases of very large energy densities, the propagating beam will essentially 
break up into smaller beams, or filaments, which severely constrains the amount of energy density that the 
beam will be able to possess as it travels through the atmosphere. Unlike the situation of passive 
propagation, the thermal blooming mechanism introduces nonlinearities into the analysis of the 
phenomena that substantially complicates a complete mathematical description. Complete analyses of 


‘See also the reference by the author, “High Energy Laser Beam Propagation in the Atmosphere: The Integral 
Invariants of the Nonlinear Parabolic Equation and the Method of Moments,’’ NASA/TM — 2012-217634, July 2012. 
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these types of propagation scenarios can only be done numerically, which was a major activity within the 
United States and Russia in the late 1980s. Other than the usual ‘order-of-magnitude’ estimates using the 
equations of fluid mechanics and wave propagation, analytical treatments appropriate for an engineering 
analysis of atmospheric propagation systems encountering thermal blooming have been lacking, making 
difficult a comparative assessment of the operation of through-the-atmosphere power transmission. 

In an effort to quickly assess the operation of power beaming methodologies proposed for application 
in continuous wave (CW) beamed energy propulsion, but at the same time, transcend the ‘order-of- 
magnitude’ analysis approach that can be employed (Ref. 3), an analytical model is needed to describe 
two very important beam parameters, i.e., beam broadening and beam steering, due to thermal blooming 
that may occur within the atmosphere in the presence of wind. (The use of the CW case represents the 
‘worst case’ scenario for power beaming since much larger power densities can be supported using pulsed 
source operation (Ref. 3).) Of course, beam broadening dictates the severity of the spread of the energy 
distribution of the beam at the target and the steering will capture an assessment of the variations of 
pointing of the beam at the target. Although such effects can be dealt with by an application of adaptive 
optics at the transmitter, an estimation of these parameters will provide a good comparative analysis of 
systems operating at different frequencies and power density levels as well as an idea of the overall 
impact of thermal blooming on system operation. (Only brief descriptions will be given concerning the 
derivation of the equations given here; the detailed derivations will be the subject of a forthcoming 
paper.) 

2.1.1 A Simple Analytical Model for Beam Broadening and Displacement 

To this end, consider the fundamental equations for a scalar electric field E(x,y,z ) propagating along 
the z-axis of a coordinate system within the atmosphere of permittivity (refractive index) s (x,y,z). 


2 ik — + V 2 E + k 2 U(x,y,z)-\)E = 0 

dx v ' 


( 1 ) 


where V]_ is the transverse Laplacian and k = 2n/'k is the wave number of the radiation of wavelength X. 

In general (Refs. 1 and 3), the permittivity s (x,y,z) is a function of the temperature T as well as a 
wavelength-dependent imaginary component s’(x,y,z) accounting for absorption and Se which accounts 
for the permittivity fluctuations due to turbulence, i.e., 

s(x, y,z) = 1 + is'(x,y,z) + sr(r) + 8e (2) 


The temperature contribution can be written as 


& t [T)-—T 
V ’ 8T 


( 3 ) 


where, for the gasses within the Earth’s atmosphere, dT/dt < 0. For the purposes of this analysis, the 
turbulent fluctuations Se will be neglected since only thermal blooming effects are being considered. (The 
combined effects of turbulence and thermal blooming should be considered in a future treatment.) Hence, 
using Equations (2) and (3) into Equation (1) gives 

r)F 

2 ik — + V]_E + k 2 (i£’(x,y,z) + e T T^E = 0 (4) 
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As discussed above, for fields possessing a large energy density, the temperature T will be a function of 
the value of the field E, i.e., T= T{E) due to the absoiption s'. This propagation equation must thus now 
be augmented by an equation from fluid dynamics that describes the creation of temperature within the 
atmosphere due to absoiption and its spatial and temporal evolution within the atmosphere. The 
appropriate equation is (Ref. 3) 


dT y2 Tl VdT _ ac 

dt 1 x dx Snxp 0 c p 

where V is the atmospheric wind velocity, taken to occur along the x-axis, a = ke' is the wavelength 
dependent absorption coefficient, c is the velocity of light, x is the atmospheric molecular diffusivity, c p is 
the specific heat at constant volume, and po is the nominal air density. (The term on the right side of 
Equation (5) is the heat source due to wave absorption and, of course, is proportional to the energy 
density \E\ 2 within the electric field.) Finally, the boundary condition giving the source of the radiation 
must be specified. In general, a good model for the field source is given by a collimated Gaussian beam 
wave at the entrance z = 0 to the atmosphere, 



Eq (x,y,0) = Eq exp 


(x 2 + y 2 ) % 


( 6 ) 


where ao is the radius of the output aperture and F is the focal length of the beam. 

Equations (4) to (6) form the basis of most investigations of the effects of thermal blooming. As 
mentioned above, the investigations can only be performed numerically for specific geometrical and 
atmospheric scenarios. In an effort to yield analytical expressions for only the most important 
performance characteristics, the following approximation procedure will be applied to these model 
equations. Assuming that a steady state condition exists for the induced temperature within the 
atmosphere, i.e., dTIdt = 0 (such an assumption cannot be made for the analysis of a pulsed laser 
situation), and taking the wind velocity to essentially move the parcel of atmosphere from the volume of 
the beam before diffusion takes effect, i.e., -VjT « (V/x){dT /dx ) , Equation (5) reduces to 



-^-\ E \ 2 

8jcp 0 Cp 


( 7 ) 


The coupled equations represented by Equations (4) and (7) will now be solved using successive 
approximations, viz, for the first approximation E\(x,v,z) for the electric field of the beam wave, 
Equation (4) will be written in which the term srT is ignored. 


dE x 

dx 


-V 2 ,E X +-E X = 0 
2k 1 2 


( 8 ) 


where the relation a = ke' which defines the absorption coefficient in terms of the imaginary part of the 
permittivity was used as noted earlier. Equation (7) becomes, within this approximation, 

dx 87TpoCp 
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Now the solution of Equation (8) is given by 


00 00 

E\ (x,y,z} - g _az / 2 J" J Eq (x',y',0)Go [x,y;x',y';z^dx'dy' 

—00 —00 


where 


G 0 (x,y;x',y',z) 


ik 

-—exp 
2 nz 


r ik 
k 2z 


({x-x'f +{y-y'f)^ 


( 10 ) 


(ii) 


is the Green function of the differential equation. Equation (8). Substituting Equations (6) and (11) into 
Equation (10) yields 


Ei(x,y,z) = -E 0 e az / 2 


f ik 

r* 2 (z)) 


Iri / q q \ 

( ik V 



exp 

— U 2 +y 2 ) 

1 + — R 2 (z) 

l z J 

V ^ J 

2 z x ’ 

V 2z )_ 


( 12 ) 


where 


with 


R2 { z ) = i \ ^f]{ l ~jR z + i Ji z ) 1 
k / 


Jr = 


1 




(13) 


(14) 


The solution of Equation (9) is written 


7) (x,y,z)- 


ac 


87ip c p V J 


J \E x (x’ ,y,z)\ dx' 


(15) 


Using Equation (12) in this relation yields, after some significant manipulation (as mentioned above, 
all the derivations that give the results displayed here will be the subject of a separate paper), 


T\ (x,y,z) = 


f \ 

a c 

SnpCpV j 


E^e~ az (A(z)) 


1/2 

f 1 ) 

1/2 

exp 

h <' V 2 1 

f 

1 + erf 

1 

1 


Uy/J 

. A (u>- 

V 

L v A ( z ) JJ 


(16) 


where 

A ( z ) = (l ~JR z f +J 2 i z2 (I 7 ) 

This ‘first approximation’ result can now be used in the second approximation of Equation (4) which 
includes Equation (16) as the thermal source, viz, 
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( 18 ) 


dE 2 

dx 


—V 2 E, +-E 2 =i-z T T\E 2 
2k 1 2 2 


the solution of which can be written 


E 2 (x,y,z) = e az / 2 J j E 0 (x',y',0)G 0 (x,y;x',y';z) 


exp 


ik 


Sj’JTi {x,y,z'}dz' 


—00 —00 


(19) 


Using this second approximation expression, one can now obtain the behavior of the intensity of the 
propagating beam wave, of arbitrary initial radius and focal length, which includes the effects of thermal 
blooming within the atmosphere. From this, one can then immediately find expressions for the beam 
broadening and beam bending. Substituting Equations (6), (11), and (16) into Equation (19) gives for the 
beam wave intensity I(x,y,L) at a point (x,y) in a transverse plane at distance L downrange from the 
transmitter 


S N 2 00 00 00 00 r 

I(x,y,L) = \E 2 (x,y,L)\ 2 =E 2 e~ az — J J J J J ex P ~y( r ( x i ' 2 -y[ 2 )- r *( x 2 -y'l 2 )) 


—00 —00 —00 —00 


+ y(( x “ *1 ) 2 + (y - ri f ) - y(( x - x 2 ) 2 + (y - y'i f ) 


, \2 \ ik 

;i) )~Yz 

+/Oi (z)(x[ -x’ 2 )-® 2 ( z ) (y\ 2 - y 2 ) 


( 20 ) 


dx\dy\dx 2 dy 2 


where 


V = y R -iyi, ®](z)s3(z) 

Tl 


T l +Y '\ ® 2 (z) = E(z)V* — ( T » +T ') 


5/2 


y 9 /2 


and 


S(z) = 


ks T E q 


aze ar l, q = ^ R 


y\ + ri 


y2p c p V j 

Again, leaving the details to a future publication. Equation (20) reduces to 


I(x,y,L) = E 2 ( x,y,L ) = El 


2 ,_,/-) 


2^-aL 


ft x (-^') fty 


exp 


(x-Av(L))“ _y2 


i L ) “H L ) 


where 

a x ( L ) = a o -^/a ( L j 

is the radius along the x-axis of the beam at a distance z from the transmitter, 


( 21 ) 


( 22 ) 


(23) 


(24) 
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a y (z) = a 0 


nl/2 


(25) 


4®^ ( z)z 2 


4 ® 2 (z)z 
k 


(l-y^z) + A(z) 


is that along the y-axis of the beam and finally, 

Ax(z) = a 0 ®i(z)z-^- (26) 

is the displacement of the initial axis of the beam along the x-direction at the target. Equation (23) can be 
used to easily calculate the energy that is transmitted to a target within the atmosphere. These relations 
indicate that in the presence of an atmospheric wind taken to be directed along the x-axis, the beam radius 
along the x-axis does not deviate from that of just pure diffractive broadening (given by Eq. (24)) but that 
along the y-axis will appreciably broaden. Elowever, the beam will also be steered along the x-axis as 
described by Equation (26). These results will now be applied to two CW power beaming scenarios to 
ascertain the comparative behavior of their overall characteristics at two disparate frequencies. 

2.1.2 Application of the Foregoing to Millimeter and Optical Wavelength Power Beaming 

Equations (24) to (26) will now be applied to assess the behavior of high power electromagnetic wave 
beams within the atmosphere at two wavelengths X of 1 .06x 1 0 6 m and 2.0x1 0 -3 m. The nominal wind 
velocity within the atmosphere is taken to be V= 4.5 m/sec (« 1 0 mi/hr) along the x-axis. For the Earth’s 
atmosphere, s t~ 2.3x10 6 K _1 , c p ~ 1.005 J (g K) -1 , and p = 1.225xl0 3 g/m 3 . The initial beam radii are 
taken to be ao = TO m. Ideally, the beams should be focused at the target, that is taken to be at a distance 
of L = l.OxlO 4 m, i.e., F = L.ln the millimeter wave case, this is realistic. That is, a millimeter wave beam 
with X = 2.0x1 0" 3 m, ao = 1.0 m at a distance of L = l.OxlO 4 m in the absence of thermal blooming will 
possess a beam radius of a(L) ~ 3.2 m (using Eqs. (17) and (24)). Elowever, at the wavelength of 
X = 1.06x10 6 m, one has a(L) ~ 3.1xl0~ 3 m (the optical wavelength beam can achieve a much sharper 
focus than its millimeter wave counterpart). At this small radii, the power density would not only surpass 
that which the theory can treat within the approximations made, but is not representative of a condition 
for optimal power transfer to the target. Thus, in the optical wave case, the focus must be backed-off, e.g., 
to achieve a beam radius at the target that is representative of that of the target size (which, for purposed 
of this preset calculation, is — 1.0 m). This dictates that for the optical case, the focal length should not be 
smaller than F = 1 .Ox 1 0 8 m. Finally, the absoiption coefficients a for each if the wavelengths must be 
found. For X = 1.06xl0 -6 m, a « 4.6x10 6 m" 1 (Ref. 4). Similarly, for X = 2.0xl0” 3 m, a « 2.2xl0" 4 m" 1 
(Ref. 5). Finally, the three power levels of Pq = 1.0x10 s W, Po = 10.0x10 s W and Po = 100x10 s W will be 
considered. (The relationship between the total beam power Po and the radiant intensity Z ( y that appears 

in Equations (22) and (23) is Pq = E^na ^ .) The operational results for such beamed transmission systems 
follows: 


Z = 2.0xl0 3 m, Z = 1.0xl0 4 m, ao = l-O m ’ F-L 
P 0 =1.0x10 6 W 

a x (Z) = 3.18 m 
a y (Z) = 4.15 m 
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A x(Z) = -1.3 m 


P 0 =10.0xl0 6 W 

a x (Z) = 3.18 m 
a y (Z) = 26.9 m 
Ax(Z) = -13.1 m 

P 0 =100xl0 6 W 


a x (Z) = 3.18 m 
a v (Z) = 267.6 m 
Ax(Z) -131.1 m 


Z = l. 06x10 6 m, Z = 1.0xl0 4 m, ao=l-Om, P = l-0xl0 8 m 

P 0 =1.0xl0 6 W 

a x (Z) = 1.0 m 
a y (Z) = 1.03 m 
Ax(Z) = -0.02 m 

P 0 =10xl0 6 W 

a x (Z) = 1.0 m 
a v (Z) = 1.34m 
Ax(Z) = -0.19 m 

P 0 =100xl0 6 

a x (Z) = 1.0 m 
a y (Z) = 4.43 m 
Ax(Z) = -1.94 m 

Clearly, the use of optical wavelengths will minimize thermal blooming effects as compared to those 
of millimeter waves. This is due simply to the fact that the absoiption coefficient is two orders of 
magnitude less at optical wavelengths than at millimeter wavelengths and, of course, diffraction effects 
dominate the larger wavelengths. The negative numbers representing the beam displacements along the x- 
axis indicate that the beam is deflected in a direction opposite to that of the atmospheric wind, i.e., the 
deflection is toward the cooler side of the beam, i.e., toward regions of larger refractive index. 
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The results of the analytical expressions derived here can also be displayed in various graphical 
formats. These will be included in the forthcoming publication that will detail the derivations of the 
equations. Equations (24) to (26) are, however, sufficient for use in the overall assessment of high-energy 
beam wave propagation within the atmosphere. 

It also remains to subject a pulsed laser scenario to the same type of straightforward analysis as given 
above for the CW case. This will form the subject of another forthcoming paper. Also, atmospheric 
conditions have been assumed in which aerosols and such do not exist that may give rise to explosive 
detonation upon heating by the radiation field. 
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